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We analyse universal statistical properties of phase shifts 
and time delays for open chaotic systems in the crossover 
regime of partly broken time-reversal invariance. In partic- 
ular, we find that the distribution of the time delay shows 
T -3/2 behavior for weakly open systems of any symmetry. 



The energy-dependent scattering phase shifts a {E) 
defined via the eigenvalues expi6 a ; a — 1, M of the 
MxM unitary scattering matrix S(E) are important and 
frequently used characteristics of the process of quan- 
tum scattering. In particular, the derivatives of phase 
shifts over energy r a = d9 a /dE are related to the du- 
ration of a collision event. For example, the quantity 
tw — M T a is the typical time delay due to scat- 

tering, the so-called Wigner-Smith time delay 111 . When 
some external parameters are taken into consideration 
(e.g. a magnetic field) the corresponding parametric vari- 
ation of the phase shifts can as well be related to some 
observables 0. 

Growing interest to the universal features Q of quan- 
tum systems whose classical counterparts demonstrate 
chaotic dynamics attracted considerable attention to the 
process of quantum chaotic scattering, see [f||j| an d ref- 
erences therein. From this point of view, different statis- 
tical characteristics of phase shifts and time delays were 
addressed in experiments on chaotic microwave reflec- 
tion H as well as in several numerical studies of various 
models of quantum scattering in disordered and chaotic 
systems 

It is interesting to mention that for the case of only two 
open channels M = 2 the phase shifts can simply be 
related to the phases of transmission and reflection coeffi- 
cients, see e.g. [^lj. The latter quantities are amenable to 
direct experimental measurements in quantum dots, see 
|l2| and references therein. Another fact attributing ad- 
ditional interest to studies of time delay statistics is that 
it is intimately connected with the issue of mesoscopic 
fluctuations of dynamic admittances of microstructures 

& 

One can extract statistical characteristics of the S- 
matrix exploiting a semiclassical periodic orbit expan- 
sion like that provided by the Gutzwiller trace formula, 
see examples of such calculations in [|]-[| ■ To proceed in 
this way, one has to employ some approximations (most 



frequently, the so called diagonal approximation). The 
resulting expressions provide an important insight into 
the problem. In particular, the semiclassical approxima- 
tion for the time delay correlations at two different ener- 
gies was derived by Eckhardt JlJ]. However, the results 
obtained in such a way have a restricted domain of ap- 
plicability; in particular, they fail to describe the system 
with only few open channels: M ~ 1. 

A powerful alternative to the semiclassical methods in 
extracting the universal (i.e. generic and system indepen- 
dent ) statistical characteristics of the scattering matrix 
is provided by the random matrix approach. In particu- 
lar, in the "Heidelberg variant" of this approach jlq] one 
relates the scattering matrix S(E) to the Hamiltonian 
of a closed counterpart of the open system. The latter 
Hamiltonian is considered to be a member of an ensem- 
ble of random matrices of appropriate global symmetry 
- an idea commonly accepted in the domain of Quan- 
tum Chaos. The use of random matrices gives one the 
possibility to apply the very powerful machinery of aver- 
aging developed by Efetov Q and to calculate different 
correlation functions explicitly for any number of open 
channels and arbitrary coupling to continua. 

When employing the Heidelberg method the actual cal- 
culation depends quite essentially on the symmetry of 
the underlying Hamiltonian. The simplest case to study 
corresponds to completely broken time reversal invari- 
ance (TRI) (systems in strong enough magnetic field), 
when the random matrix Hamiltonian H is taken from 
the Gaussian Unitary Ensemble (GUE). For such systems 
the statistics of phase shifts, delay times and resonance 
poles was thoroughly investigated by two of us recently 
p7[ , see a detailed exposition of the calculation in jjl8|. 
In the opposite case of fully preserved TRI when H is 
a member of the Gaussian Orthogonal Ensemble (GOE) 
some aspects of time evolution of a chaotic system were 
considered in |2(|, the correlation function of Wigner- 
Smith time delays for two different energies was found in 
]l9t and the distribution of time delays was obtained for 
the perfect coupling case in Jl3| , see also . Let us also 
mention the paper |22j addressing the issue of parametric 
correlations for S— matrix elements. 

In the present communication we extend the analysis 
of statistical properties of phase shifts and time delays 
to the whole crossover region of gradual breaking of the 
TRI. Different characteristics of chaotic and disordered 
systems in this crossover regime were under quite an in- 
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tensive theoretical investigation recently 23 - 25|] . 

Within the framework of random matrix theory, 
Hamiltonians of the closed chaotic systems under consid- 
eration are conveniently represented as |^,^|: H(y) — 
H s +i-JL=H A , where H s is N x N GOE matrix and H A is 
a real random antisymmetric matrix of the same dimen- 
sion. For the sake of generality the symmetric matrix Hs 
is taken in the form @ : H s = H { s ° ] + -^H^ . This form 
allows one to simulate the influence of such perturbations 
(e.g. a variation of the strength of scattering potential) 
which do not break the TRI. All elements of random ma- 
trices are independent and normalised in such a way that 
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= N. 
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In the limit N — ► oo the crossover is driven by the 
parameter y 6 [0, oo), with y = oo corresponding to 
completely broken TRI. Physically the parameter y is 
proportional to the magnetix flux through the system <E>. 
One may also notice that the typical shift of the levels 
due to the antisymmetric perturbation is SE y /A ~ y 2 
pEfl , where A is the mean level spacing. 

Within the framework of the Heidelberg approach 
Jl5[ the coupling of the chaotic region to the incom- 
ing/outgoing waves is described with the help of the M x 
N matrix W of amplitudes Wi a , a — 1,2,...,M; i — 
1,...,N, which couple the internal motion to M open 
channels. In what follows we consider the case of arbi- 
trary, but fixed M whereas N — > oo. Without much loss 
of generality these amplitudes can be chosen in a way 
ensuring that the average S— matrix is diagonal in the 
channel basis: {S a b} = o~ab{S a a}- The strength of cou- 
pling to continua is convenient to be characterized via 
the " sticking probabilities" (also called the " transmission 
coefficients") T a = 1 — KS^a)! 2 which are given for the 
present model by the following expression Jl5| ]: 
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with v(E) = 7r _1 (l — E 2 /4) 1 / 2 being the density of states 
for the GOE matrices related to the local mean level spac- 
ing as A = (vN)^ 1 . The quantity T a measures the part 
of the flux in channel a that spends a substantial part of 
time in the interaction region jL5| . We see that both lim- 
its 7 a — > and 7 a — > oo equally correspond to the weak 
effective coupling regime T a « 1 whereas the strongest 
coupling (at fixed energy E ) corresponds to the value 
7 a = 1 . The maximal possible coupling corresponding to 
the upper bound T a — 1 is achieved in the present model 
for an energy interval in the vicinity of the center E = 0. 
Below we restrict our attention to this point of spectrum 
in order to present our final results in the most compact 
form. Moreover, we consider all channels to be statisti- 
cally equivalent: i a — 1 for a = 1, M. Generalization 
to arbitrary E and non-equivalent channels can easily be 
done, see |l8[ . 

In the earlier work fl7],[l^] it was shown that one can 



study very effectively the statistics of phase shifts 9 a con- 
sidered modulo 2ir. To this end we find it to be conve- 
nient to introduce the auxiliary " phases" <j> related to the 
phase shifts 6 as <fi — arctan{7 _1 tan (0/2)} and consider 

the density pE,x,y{4>) = ^Ea^W"^^. 1 ^))- 

The connected part of the correlation function of these 
densities is our main object of interest. It can be found 
performing calculations mutatis mutandis, similar to that 
presented in p^JT^ ] and it turns out to be dependent only 
on the difference <f> = <f>i — </> 2 : 
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(</>) = (pE=0,x=0, V i(<l>l)PE=n,x,yi(<h))c 



ReJdxJdMj^M^e-^ 2 ^-^-" 2 ^ x 
-li i 

e - l u 1 (X 1 X 2 -\)+y 1 y 2 (Xl-Xl)-^{y 2 1+ yl){\l + \l-\ 2 -l) y 

{(1 - A 2 ) cosh a - (A 2 - A 2 )sinha+ 
ft [ivl + yl) (C 1 - A 2 ) cosh a + (A 2 . - A 2 )sinha) + 
2yiy 2 (A 2 + A 2 + A 2 -l)sinha]} (2) 

where the "channel factor" is equal to 
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1 + 2iAi A 2 tan 4> - tan 2 4>(X\ + A§ - 1) 
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(,g + A 1 A 2 ) 2 -(A 2 -l)(A 2 -l) 



here g — 2/T — 1, with T being the transmission coeffi- 
cient introduced above. 

Secondly, one can extract explicitly the general distri- 
bution function of the scaled partial delay times V(t s ) = 
(if J2a S ( T s ~ ^ T a(y))) in tnc crossover regime: 
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where C M = [{2n) l / 2 2 M / 2+1 T (f + 1)] and 



Ji(A 2 ) = JodM^^e-^hM^ 
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J 2 (X, A 2 ) = Ay 2 [(1 - X 2 )e-P + A 2 (l - e^)] - (l - e 
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and where lZ,a,u) denote 1Z = A 2 
1, a — 2/11/2(1 — A 2 ) and u> = nfl/A. 

The correlation function presented above is a very in- 
formative object. First of all, having it at our disposal 
it is a relatively easy task to show [Q that the correla- 
tion function of Wigner-Smith time delays: K^ x yi y2 — 
(t w (E = 0,x = 0, yi )T W (E + n,x,y 2 ))/(T W {E)) 2 is 
given in the crossover regime by the same expression 
Eq. (||) provided one replaces the " channel factor" Eq. (j|) 
by: 



M/2 
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where v(ip) = g — y/ g 2 — 1 cos ip, (3 = 2y 2 (l — A 2 ) and 
I p (z) stands for the modified Bessel function. 

The distribution Eq. (||) is valid for any number of open 
channels M and any value of transmission coefficient T 
and as such is quite complicated. To get a better un- 
derstanding of its typical features it is reasonable to look 
separately at two limiting cases of strong/ weak coupling 
to continua. 

For the strong coupling regime T — 1 (i.e. g = 1) 
the ip— integration in Eq.(||) drops out, but the result- 
ing expression is still quite cumbersome. However, one 
can easily find the long time asymptotics to be of the 
following form: 

V(r » 1) - / U M {y)r7 (2+M) y>0 () 

where Um{u) is a rather complicated function of the sym- 
metry breaking parameter y. 

The second of this expressions holding for unbroken 
TRI is actually exact for arbitrary r s as can be seen per- 
forming the limit y —> in the general Eq.(|5|) at T = 1. 
This fact was first conjectured in |2l]] and derived for the 
particular case M = 1 by another method in |Tg| ]. 

Eq.(||) demonstrates that the limits y — > and r — > oo 
do not commute. To understand this phenomenon better 
it is instructive to consider the case of " weakly broken" 
TRI, y <C 1, in more detail. A close inspection shows that 
for such regime there emerges one more relevant time- 



scale 



'yocy 



2 ^> 1 such that for the domain 1 <C r s -C t v 



whereas at t s 3> r y this be- 

-(M+2) 



the distribution function V(t s ) shows the GOE-like be- 
haviour: V(t s ) CX T s (- A/ / 2+2 ^ 

haviour changes to the GUE-like: V(t s ) cx T y ^ 2 , s 

We suggest the following transparent physical interpre- 
tation of the scale t v : this is just the time H/5E y nec- 
essary by the Heisenberg uncertainty relation to resolve 
a typical shift 5E y due to the TRI-breaking perturba- 
tion. If the particle dwells in the scattering domain for a 
time shorter than r y it can not "feel" the magnetic field 
effects and the corresponding asymptotics is GOE-like. 
However, for large enough times the particle resolves the 
effect of broken TRI whatever small is the magnetic field. 
This explains why the most distant asymptotics of the 
time delay distribution is always GUE-like, provided the 
magnetic field is not identically zero. 

Let us now turn our attention to the opposite limit of 
an almost "closed" chaotic system: T « 1. Exploiting 
j> lwe find for arbitrary number of open channels and 
arbitrary y the following universal (up to a coefficient) 
behaviour of delay time distribution: 



V{t s ) (x g~ 1/2 T~ 3/2 when 



g L <.T a 4Zg (7) 



in the parametrically large region of delay times. The 
proportionality coefficient in this formula depends on the 
symmetry breaking parameter y and M in a complicated 
way. 



Such a t~ 3 / 2 behaviour holding irrespective of the TRI 
symmetry is the most robust feature of the time delay 
statistics of weakly open chaotic systems. It was first 
obtained in |lq] for the case of broken TRS, but phys- 
ical arguments show that it is a very general feature 
simply following from the picture of well-isolated reso- 
nances typical for such systems [Q . It also seems to be 
quite insensitive to the particular details of definitions of 
time delays and holds for distributions of such slightly 
different quantities as Wigner-Smith time delay, partial 
delay times or even "dwell times". In Fig.l we plotted 
a typical fluctuating pattern of energy-dependent "dwell 
times" obtained in the paper JTo| ] in the course of numer- 
ical simulations of quantum chaotic scattering on a two 
dimensional cavity in tunneling contact with two waveg- 
uides. Sampling the distribution of dwell times over the 
chosen range of energies we find a good agreement with 
the predicted t -3 / 2 behaviour, see Fig.l. 




FIG. 1. The distribution of the dwell times r in chaotic 
scattering in weakly open Sinai-like billiard. Inset shows a 
fluctuating pattern of dwell time versus energy (data were 
kindly provided to us by the authors of Ref. [10]). Solid line 
shows the theoretical prediction V(t) <x t -3//2 . 

Outside the parametrically large interval g~ x <C r s <C 
g our general expression Eq.(^|) predicts an exponential 
cutoff at t s < 9 and a crossover to the behaviour 
described by Eq.(|^) for the asymptotically large times 
r 3> g. One can check that for "weakly broken" TRI 
again there emerges a scale r y (g) ~ g/y 2 such that the 
asymptotic tail is GOE-like at g <C r s <C T y (g), but al- 
ways GUE like for r s 3> r y (g) , in full agreement with the 
discussion above. 

Finally having in mind the comparison with the semi- 
classics let us consider in more detail the large-channel 
limit Jlf > 1 of our general expressions like Eq.(||) 
describing the correlations of phase shifts and Wigner- 
Smith time delays. When doing this it is natural to con- 
sider the angle difference (j> = <f>\ — <p2 to be of the order 
of 4> ~ 1/M <C 1. Then one rescales = <p/M, substi- 
tutes tan0 ~ (f>/M in the "channel factor" Eq.(||) and 
performs the limit M — > oo explicitly. For lu = the re- 
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suiting expression turns out to be identical to the para- 
metric correlation function of eigenvalues of large ran- 
dom matrices in the crossover regime derived for the first 
time by N.Taniguchi et al. f25|| . Taking into account that 
<t>\ — (f>2 ~ 1 /M results also in 61 — 62 ~ 1/M, we conclude 
that the statistics of scattering phase shifts in the large 
M limit is just the same as that of energy levels of closed 
chaotic systems. The latter conclusion is in agreement 
with the available numerical results obtained for a real- 
istic model of chaotic systems with M = 23 in |J. It is 
also interesting, that the only modification required for 
lu is to replace <p — ► cf> + lo. 

Considering the time delay characteristics one should 
take into account that the width of the time delay distri- 
bution is of lower order in M as compared with the mean 
value (tw) when resonances are overlapping: MT 1 
fl9f . To extract the time delay correlations in the corre- 
sponding limit requires a calculation similar to that done 
in the paper by Pluhar et al. (2^]. The resulting expres- 
sion turns out to be quite a transparent one and is given 
by: 
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T i - LU 



[r? 



(8) 



provided that T± = MT/2 + x 2 + {y x ± y 2 ) 2 > 1. Ac- 
tually, this formula is nothing else but the semiclassical 
expression for the time delay correlator. It can be ob- 
tained from the Gutzwiller trace formula in diagonal ap- 
proximation, with the quantity MT/2 being replaced by 
the decay rate out of the chaotic region, see JlJ|,[l8| • 

The last point to be mentioned is related to the issue 
of fluctuations of low-frequency admittance as defined 
in We noted above that the time delay fluctuates 

weakly in the many-channel limit M — > 00. Using this 
fact and the relation between the time delay and the 
low-frequency admittance G 1 (u>) presented in [|l3| (see, 
however, |p6| ) one finds that the parametric correlator 
(G I (0 1 yi)G~'(x,y2))/(G 1 ) 2 - 1 of the admittance in the 
limit M 3> 1 is given by: 
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(9) 



where C e denotes the so called "geometric capacitance" 
p3[ and e stands for the electron charge. In the limiting 
cases of unbroken (yi — y2 — 0) and completely bro- 
ken (yi = t/2 - * 00) TRI this expression coincides with 
that found recently by another method by Brouwer and 
Buttiker |J|. 
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